Rectification of spatial disorder 
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We demonstrate that a large ensemble of noiseless globally coupled-pinned oscillators is capable of 
rectifying spatial disorder with spontaneous current activated through a dynamical phase transition 
mechanism, either of first or second order, depending on the profile of the pinning potential. In the 
presence of an external weak drive, the same collective mechanism can result in an absolute negative 
mobility, which, though not immediately related to symmetry breaking, is most prominent at the 
phase transition. 
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Rectifiers are special devices capable of extracting a 
steady output signal (current) even from a perfectly 
center-symmetric input signal (drive) [lj]. Their min- 
imal operating conditions typically require an asym- 
metric internal dynamics, possibly chosen to optimize 
performance, and a non-stationary unbiased input sig- 
nal, mostly a cyclostationary periodic drive or a time- 
correlated noise. In thermodynamical terms, rectifiers 
are devices that operate under nonequilibrium conditions. 
Direction and magnitude of their output current strongly 
depend on both the drive(s) (intensity and time scales) 
and their intrinsic noise (temperature). Rectifiers are of- 
ten called ratchets or else, to emphasize the role played 
by fluctuations, Brownian motors 0]. 

Rectification of an external symmetric signal can be 
strong enough to overcome the action of an additional 
static drive, or load, so that the rectifier can do mechani- 
cal work against it. Moreover, under appropriate circum- 
stances, the load itself induces an asymmetry in the oth- 
erwise perfectly symmetric dynamics of the motor, thus 
allowing rectification of symmetric ambient fluctuations 
in a direction that depends on the details of the modified 
motor dynamics. This mechanism can produce the ap- 
parently counterintuitive effect of a spontaneous system 
response oriented against the load: Rectifiers working un- 
der such operating conditions are said to exhibit absolute 
negative mobility (ANM) @|. 

The question now arises whether only time-dependent 
signals can fuel such rectification mechanisms. To ex- 
plore the possibility of replacing temporal variability by 
spatial randomness, we considered an ensemble of N cou- 
pled simple phase oscillators, each driven by a static 
force, randomly chosen to model their interaction with a 
spatially homogeneous, quenched disordered landscape. 
Their equilibrium phase values coincide with the local 
minima of a periodic bistable pinning potential with 
center-symmetric cells. Following Kuramoto's prescrip- 
tion strong cooperative effects among oscillators are 
ensured by assuming global pair coupling, irrespective 
of spatial separation. At finite temperatures and in the 
limit N — !• oo, this system is known to undergo sponta- 



FIG. 1: (Color online) Phase diagram at f e — 0. The thresh- 
olds of SSB (K c ), global bistability (K w ), and ANM (K m ), are 
plotted versus b when a u = a = 1. All curves are obtained by 
analyzing Eqs. ([3])-([6| in the mean field approximation. Nu- 
merical data (dots) are obtained by integrating Eq. |TJ with 
N — 10 6 , starting from random initial conditions. Note that 
the two transition branches (Kf. and K^ 1 ) are connected by 
a cusp at (K*,b*) ~ (1.68,0.94). The dashed curves are the 
analytic limits of K and K^ 1 for a u — 0, i.e., g(uj) = 5(uj)\ 
the horizontal arrow indicates the critical coupling for the 
unpinned Kuramoto model, ^8/%^. 



neous symmetry breaking (SSB) Q. As a consequence, 
it generates a spontaneous current in either direction; 
moreover, subjected to an external weak global drive, it 
can also exhibit ANM and thus work as a motor. Since 
first detected, both effects have been interpreted as rec- 
tification phenomena induced by thermal noise 0, 0] • 

In this Letter, we show that a large ensemble of noise- 
less coupled-pinned oscillators is capable of rectifying lo- 
cal forces or stresses, associated with spatially distributed 
quenched disorder. By numerically analyzing the station- 
ary dynamics of the oscillator ensemble, we demonstrate 
that the onset of spontaneous rectification currents obeys 
a dynamical phase transition mechanism, either of first 
or second order, depending on the degree of bistability of 
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the pinning potential. Moreover, we observe that ANM 
is neither induced by thermal noise, nor immediately re- 
lated to symmetry breaking, although its magnitude is 
most prominent at the phase transition. 

We consider a system of N pinned phase oscilla- 
tors coupled through the Kuramoto-type global interac- 
tion 

N 

& = Ui + f e - v'ifa) - (k/n) - fa), (i) 

where <pi and Wi denote, respectively, the phase and the 
intrinsic frequency of the «-th oscillator. The intrinsic 
frequencies are assumed to be randomly distributed ac- 
cording to a Gaussian distribution function, g(ui), with 
zero mean and variance <r u ] a tunable frequency bias is 
introduced by the external drive f e . The third and fourth 
term on the r.h.s. represent the pinning force acting on 
the i-th oscillator and all-to-all ferromagnetic (K > 0) 
coupling between oscillators, respectively. Due to the 
global nature of the interactions, mean-field (MF)-type 
phase transitions are expected. 

For simplicity, the pinning potential V{<f) is chosen to 
be identical for all oscillators, 



V(cj)) = -a cos (f> + (6/2) cos 2<j>, 



(2) 



with V'(4>) = dV/d<f>. V{4>) is a symmetric biharmonic 
potential, V{—<j>) = V{4>), with period 2tt. Its unit cells 
are monostable with one minimum at (j> = (mod 2tt) 
for b < a j 2, and bistable with two stable points at <j) = 
±<f) m for b > a/2, where <f m = arccos(a/26). Isolated 
oscillators are symmetrically locked for + f e \ < f p 
with locking phases |</>j| < (f> p , where the largest locking 
phase 4>p is the (larger) positive solution of the equation 
V"{4>) = and the depinning frequency threshold is f p = 

V(4, p ). 

In contrast to Refs. [y, [Z(, the coordinates <f>i are not 
coupled to a heat bath. As the ergodic assumption does 
not apply here, the oscillator dynamics can depend on 
the initial conditions (i.e.), {^(0)}. The simulation data 
reported below refer to the case of disordered i.e., where 
{(f>i{0)} have been uniformly randomized in the interval 
(0, 27r). We also simulated ordered initial configurations 
with pinned phases, 4>i(0) = m , and even annealing pro- 
cedures. We found that the key qualitative features of our 
main conclusions do not change. 

Following the Kuramoto's approach we introduce 
the synchronization order parameter A and the aver- 
age phase 6>, defined by Ae t6 = (e 4< ^), where (Oj) = 
(1/-AQ X^j=i Oj denotes the oscillator ensemble average. 
Accordingly, Eqs. ([lj can be rewritten as 



^ = w i +/ e -[/ / (^), 
with the effective potential U {(f) given by 
U{4>) = V{4>) - KAcos{<j) - 6] 



(3) 



(4) 
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FIG. 2: (Color online) Spontaneous currents. \J\ (squares), A 
(dashed) and |5*| (solid) versus K (a) in region I at b = 0.75, 
and (b) in region II at b = 1.5, from the numerical integration 
of Eq. ©. 



Here, the parameters A and are determined through 
the the self-consistency relations 

C = A cos 8 = {cos 4>j) = J cos {cj){u})) g{uj)duj, (5) 

S = A sin 8 = (sin 4> 3 ■) — I sin {<fi{uj)) g{uj)du). (6) 



For 9 = 0, the potential U{(f>) is characterized by the 
depinning frequency thresholds ±/ p , the corresponding 
depinning phases ±</> p , and, if bistable, the two symmet- 
ric minima at (j) = ±0 m , similar to the isolated oscil- 
lators. The interaction term in U{(f) tends to suppress 
the bistability of the pinning potential even for b > a/2: 
U{(f>) is bistable only for small K, < K < K w , where 
the bistability threshold K w is the solution of the implicit 
equation of K W A{K W ) = 2b — a. For 8^0, U{<p) is no 
longer mirror symmetric; correspondingly, the symmetric 
depinning thresholds, ±/ p , are replaced by two distinct 
thresholds, f±{KJ e ), with /+ ^ -/~. 

In the absence of external drives {f e = 0), the anti- 
symmetry of Eq. (|3|) guarantees 4>{uj) = — (/>(— uj) as long 
as 8 = 0, so that the current J = {(j>)/2ir is identically 
zero. A net SSB current sets on only for 8^0, where two 
distinct frequency thresholds {f p ^ — fp) break the left- 
right symmetry of the running oscillators. As the pinning 
ensures a certain degree of synchronization (A > at all 
K > 0), the average phase 8 alone (or more conveniently 
S) becomes the proper order parameter of the system. 
Supersymmetry considerations [l[ rule out spontaneous 
currents for a purely harmonic pinning potential (a = 
or b = 0). A biharmonic V{<fi), instead, allows nonzero 
S and J in the strong coupling regime of K > K c . The 
transition threshold K c is actually a function of two pa- 
rameters only, a u and 6, since a can be rescaled to unity 
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without loss of generality. 

To investigate a large ensemble of disordered oscilla- 
tors, we had recourse to the numerical integration of the 
equations of motion, Eq. ([T]), for the entire ensemble. For 
the sake of a comparison, we also numerically solved the 
self-consistency equations, Eqs. ([3])- ([6]). In Fig. []] we 
compare the dependence of K c and K w on the pinning 
bistability parameter b at f e = 0. The intersection of 
K c (b) and K w (b) at b = b* defines two distinct dynami- 
cal regimes: (i) region I (a/2 < b < b*). The SSB tran- 
sition occurs after the pinning bistability is completely 
suppressed by the Kuramoto coupling (K^ > K w ). The 
transition branch K c diverges as b — > a/2 as expected 
and decays toward a minimum at b = 6*; (ii) region II 
(b > b*). The transition branch K^ 1 always lies below 
K w [1] , so SSB occurs inside the bistability regime where 
multiple solutions are possible. Moreover, we notice that 
K*(b) and K^ 1 ^) form a cusp at b = b* , which indicates 
that two different SSB mechanisms are at work in regions 
I and II. 

Following this lead, we numerically computed the 
spontaneous current, J, as a function of K at b = 0.75 
in region I and b = 1.5 in region II, see Fig. [21 The dif- 
ference is revealing: The onset of J can be regarded as a 
dynamical phase transition of the second order, with 

\J\^(K-K c f'\ (7) 

in region I, and of the first order, in region II. In the 
latter, |J| jumps discontinuously from to a maximum 
at K = K* 1 , and finally decays exponentially at larger 
K. Similar continuous or discontinuous behaviors at the 
transition points are exhibited by the corresponding or- 
der parameters A and S. No first order transitions were 
detected in the presence of noise . 

The existence of a nonzero SSB current for K > 
can be analytically explained in region I by linearizing 
Eqs. ([I])-© for small S ~ A9. A net current can only re- 
sult from the oscillators running respectively to the right 
with ui > /+ and to the left with oj < f~ . In the linear 
regime, it is easy to find from Eqs. © and (UJ) 

f±~±f p + K\cos4> p \S, 

which yields the net current J for small drive f e 

2irJ oc -(/+ + f~)/2 ~f e -K\ cos&IS. (8) 

Here, Eq. ([5]) implies a linear response also to the external 
drive / e , which will be discussed later. At f e = 0, J oc — S 
near the transition. 

The critical behavior of S near the transition can be 
also extracted by linearizing the self-consistency Eqs. (O 
and © at f e = 0. Since the system has the Ising-like Zi 
symmetry, Eq. ([6]) is expanded in odd powers of S as 

S = {K/K C )S - kS 3 + 0(S 5 ). (9) 



With k > 0, we find the stable solution of S = for K < 
K[ = K c and a new stable branch of S oc (K - K') 1 / 2 
for K > K^,, which determines the MF critical exponent 
of the current in Eq. ([7]). 

The first order transition along K^, 1 (b) in region II 
is a unique feature of our model. It occurs where the 
static phase solution of Eq. ((3|) for the locked oscilla- 
tors consists of two disconnected branches, 4>i,2 (w), with 
<pi (-co) = — 02 (^) in a small lo interval centered around 
u) = 0. For a homogeneous randomization of the initial 
phases, both solutions </>i,2(w) contribute to the ensem- 
ble averages with statistical weights proportional to the 
respective basin size. For oscillators with u) = 0, these 
weights are equal at 9 = 0. As 9 departs from zero, the 
symmetry of two branches is broken, but the determinis- 
tic nature of dynamics does not allow the system to redis- 
tribute the oscillator ensemble across the gap separating 
two solutions (as it would in the presence of noise @). 
Such a resistance of the system against perturbations re- 
flects itself in a delayed onset of the SSB phase in region 
II; as a consequence, k becomes negative well before the 
transition, which causes a discontinuous jump in 9 at the 
transition. Indeed, the SSB transition is delayed until the 
5 = solution becomes locally unstable (no matter what 
the global potential shape is, as the oscillators are un- 
able to redistribute themselves between the two basins). 
Therefore, the transition threshold K^, 1 is again deter- 
mined by the linear term in Eq. ([9]), with the difference 
that for k < (at the transition) the stable solution of S 
exhibits a discontinuity. This mechanism is distinct from 
the ordinary first order transition in equilibrium systems. 

In the regime of weak disorder, the con- 

tributions from the running (depinned) oscillators to C 
and 5* are negligible. Then, the ensemble averages in 
Eqs. (J5J) and ^ can be restricted to static (pinned) os- 
cillators only with w € (fp,fp)- In this approxima- 
tion, we can locate K c in Eq. ([9]) very accurately, leading 
to the curves K[(b) and K^lb) drawn in Fig. [TJ One 
may obtain a simpler but less accurate expression for 
K c as 1/K C ~ (cos 2 (/)/U"(4>)) Q, which becomes ex- 
act at a u = for identical oscillators Its solutions 
exhibit, besides the two diverging branches for b — » oo 
and b — > a/2, also a suggestive lower bound at b = b* 
as 0-^^%/tt, which happens to coincide with the critical 
coupling of the unpinned Kuramoto model [V(4>) = 0]. 

The reentrant profile of the transition curve K c (b) can 
be explained through a simple qualitative argument. The 
dynamics of a pinned oscillator in Eq. (TTJ results from the 
competition of two opposing forces: (i) the random lo- 
cal drive with typical amplitude a u ] (ii) the Kuramoto 
coupling, which is enhanced by a bistable pinning po- 
tential. For a phase transition to occur, the average 
attractive Kuramoto force must win over disorder, i.e., 
(K/2)$m2(j) m = K{a/2b)y/l - (a/26) 2 must be larger 
than o~ u . The non-monotonicity of sin2</> m in the range 
< <\>m _ t/2 determines the convexity of K c (b). In par- 
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FIG. 3: (Color online) Absolute negative mobility, (a) J-f e 
characteristic curve for b = 0.75 and three different K from 
numerical integration of Eq. (JTJ. (b), (c) fio versus K in 
the vicinity of the transition point for (b) b — 0.75, and (c) 
b — 1.5. For graphical convenience, (Mo has been rescaled by 
the factor 2n/{\<f>\). 



ticular, K c {b) diverges in correspondence with the zeros 
of sin2</> m , like |(& — a/2) -1 / 2 for b — > a/2, and propor- 
tional to 2b for b — > oo, in agreement with Fig. [1] 

Further evidence of the coexisting phase transitions of 
the first and second order was obtained by investigating 
the system response to a finite bias f e , and more specifi- 
cally, by analyzing the driven current J(/ e ) and its zero- 
point mobility, fio = {dJ/df e )f e =Q. We computed both 
quantities by direct integration of Eq. (JT]) and summa- 
rized our results in Fig. |3l In panel (a) , the characteris- 
tic curve J-f e clearly exhibits three different regimes: (1) 
< K < K m . Below a certain threshold K m (also plot- 
ted versus b in Fig. [TJ , J is parallel to f e as expected in 
the linear response theory ((j,q > 0); (2) K m < K < K c . 
The drive modifies the global interaction of the oscillators 
with their pinning potential in such a fashion that the 
ensemble response points against the drive. The curve 
J(/ e ) is continuous at f e = with (j,q < (ANM), and 
turns positive for larger / e ; (3) K > K c . The slope of 
J(/ e ) at the origin, (iq, can grow so negative that eventu- 
ally the curve splits into two disconnected antisymmetric 
branches with J(0 + ) = — J(0~) < and uniquely defined 
negative (1q (ANM). This discontinuity is a signature of 
SSB and is a common feature of both regions I and II. 
Note that ANM occurs regardless of the presence of ther- 
mal noise. 

The different SSB mechanisms of regions I and II also 
influence the ANM properties of the system, see Figs. [3] 
(b) and (c). In region I, (J.q(K) develops an asymmet- 
ric negative peak numerically compatible with a two- 
sided divergence for K — > K^. In region II, (io{K) 
is clearly discontinuous with a negative diverging branch 
for K — > Kj, 1 ^, and a slowly decaying one for K > K]} . 



In both regions, for exceedingly large or small coupling 
constants (io{K) tends to vanish, as expected. Finally, we 
remark that ANM does not necessarily anticipate SSB, 
as apparent in Fig. [IJ where the curve K m (b) (marking 
the ANM onset) crosses into a region of the monostablc 
pinning for b < a/2, inaccessible to K c (b). 

The numerical results of Fig. [3] can be qualitatively un- 
derstood from the linear response Eq. ([8]) for the current, 
which yields 

ex 1 - K\ cos^p|x, (10) 

where \ — (dS/df e )f e —o is the zero-field susceptibility. 
The curve K m {b) plotted in Fig.[T]was obtained, indeed, 
from Eq. (fTUl) and then checked against the numerical 
integration data from Eq. ((T|). Note that K m (b) also 
develops a minimum (actually a cusp) as it crosses K w (b) . 
Like in the ordinary equilibrium MF theory, one can find 
X oc \K C — K\^ 7 (7 = 1) near K = K\ and, notably, also 
for K < K^ 1 , where the local instability occurs around 
5 = 0. Therefore, ANM also diverges at SSB transition 
point with MF susceptibility exponent 7 = 1, which is 
consistent with our numerical results. 

Finally, among the possible applications of our study, 
we mention the tug- of -war between two groups of molecu- 
lar motors [ll|. Indeed, molecular motors running to the 
right and to the left, due to a Kuramoto type of attrac- 
tion, can pull polarized filaments in opposite directions. 
However, to treat mobile motors in space, the interaction 
topology of our model needs to be modified ad hoc [IH • 
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